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Abstract. The spectroscopic data for the short-period (0!l984) eclipsing binary V38, discovered by the OGLE micro-lensing 
team in Baade's Window field BW3, are analyzed. Radial velocity curves are derived from mid-resolution spectra obtained 
with EMMI-NTT at ESO - La Silla, and a simultaneous solution of the existing light curve by OGLE and of the new radial 
velocity curves is obtained. The system is formed by almost twin M3e dwarf components that are very close, but not yet in 
contact. The spectra of both dwarfs show signatures of the presence of strong chromospheres. Spectroscopy definitely confirms, 
therefore, what was suggested on the basis of photometry: BW3 V38 is indeed a unique system, as no other similar binary with 
M components and in such a tight orbit is known. 

Within the limits posed by the relatively large errors, due to the combined effect of system faintness and of the constraints on 
exposure time, the derived physical parameters seem to agree with the relations obtained from the other few known eclipsing 
binaries with late type components (which indicate a discrepancy between the available evolutionary models and the data at 
~ 10% level). A possible explanation is the presence of strong magnetic fields and fast rotation (that applies to the BW3 V38 
case as well). A simple computation of the system secular evolution by angular momentum loss and spin-orbit synchronization 
shows that the evolution of a system with M dwarfs components is rather slow, and indicates as well a possible reason why 
systems similar to BW3 V38 are so rare. 
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,£h 1- Introduction 



OGLE BW03 V038 ( BW3 V38, for brevity, through this pa- 
per; a 20 oo = 18 h 04 ra 44 s :i9, <5 2 ooo = -30°09'05", P = 0?19839, 
I max = 15.83, (V - I) max = 2.45, AI = 0.78) appears in the 
second in stalment of the OG LE Catalog of periodically vari- 
able stars (Udalski et al. 1995) and is therein classified, on the 
basis of light curve shape, as a contact binary of W UMa-type 
(EW). It was later rejected, however, from the bona-fide sample 
of contac t syst ems in the OGLE fields, which was selected by 
iRucinskil dl997h by means of an objective Fourier-based clas- 
sification of the light curves. The system was found marginally 
outside the region that, in the space of parameters defining the 
light curve properties, corresponds to the contact configuration. 

Maceroni & Rucinski (1997, hereafter paper I) analyzed 
the OGLE photometric data, providing a solution of the I-band 
light curve and studying in detail the related uncertainties. They 
reported, as well, that the system is undetected in the ROSAT 
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survey, with an upper limit to the X-ray flux of 1.8 x 10~ 13 
erg crrT 2 s _1 (this result was not surprising on the basis of the 
system properties and of a rough estimate of its distance of 
~ 400 pc). 

The de-reddened V-I color of BW3 V38 ((V - I) = 2.3 - 
2.4, paper I), its light curve solution and the short period sug- 
gest a configuration with two similar very close and fast rotat- 
ing M-dwarfs. This makes of BW3 V38 an extremely interest- 
ing object on several grounds. 

While a few eclipsing binaries with M-t ype componen ts are 
known at present (see the recent paper by Ribas (2003) for a 
summary) BW3 V38 is the only one with so tight an orbit, and 
this, according to the current knowledge on close binary forma- 
tion, is very unlikely to be its original configuration, as "in situ" 
formation is ruled out for such small separations (see Section 
4). 

We think, therefore, that BW3 V38 started its MS life with 
a larger separation and has evolved through the combined ef- 
fect of angular momentum loss by magnetic braking (AML) 
and of spin-orbit coupling. The presence of a strong magnetic 
field was suggested indeed by the typical magnetic activity sig- 



2 



C. Maceroni and J. Montalban: Absolute elements of BW3 V38 



natures found in the photometric data, as the solution of the 
OGLE light curve (paper I) required a large dark spot on the 
stellar surface of the primary component. We expected, there- 
fore, to detect other typical activity indicators in the spectra. 

In synchronized binaries the net effects of AML are shrink- 
age of the orbit, orbital period decrease and faster spinning of 
the components. The secular evolution of the orbital period is 
strictly related to the efficiency of magnetic braking. 

At the time when the available data were rather sparse and 
inhomogeneous, the lack of close binaries with late (M) type 
components could be explained in terms of a selection effect 
against the discovery of intrinsically faint objects. However, 
the modern large surveys for variability (such as OGLE) have 
shown that the effect is real, as these systems are indeed very 
rare, even in large, homogeneous and at some level complete 
datasets. 

The dearth of such systems implies then some evolution- 
ary effect. Either, for some configuration, the orbital evolution 
should become so fast to prevent system detection or, on the 
contrary, it should be very slow. The latter hypothesis seems 
more likely in light of the current understanding of the AML 
processes in late-type s ingle stars. 

Studies of cluster jTerndrup et alJ l2 000l and references 
therein) and of field M dwarfs ( Delfosseet al. I ll 9981) of differ- 
ent ages favor the hypothesis of a spin-down time-scale that 
increases with decreasing mass (and show a saturation of the 
activity indicators with rotation). This evidence has been re- 
lated to the properties of the aQ. dynamo that operates at the 
interface between the radiative core and the convective enve- 
lope, and that is expected to work up t o spectral type M3 (or 
maybe later for magnetic dwarfs, see Mullan & MacDonald 
The dynamo and the braking efficiency increase - 
at fixed radius - with faster rotation, until a "saturation" of 
the mean field is reached at rotational velocities that slightly 
depend on the sample, but are of the order of 10 km s" 1 . 
Looking instead at the behavior with spectral type, as the angu- 
lar momentum is essentially extracted from the convective en- 
velope, which deepens with advancing spectral type, later stars 
have longer spin-down timescales. Finally, in fully convective 
stars a further reduction of AML is expected, as the dynamo 
ch anges in nature (perhaps the small-sca le turbulent dynamo 
of lDurnev. De Young & Roxburg h! dl993h ^ . 

In synchronized binaries with late type components, that - 
at variance with single stars - rotate faster with increasing age, 
both effects contribute to lengthen the "contact" timescale. The 
studies of the period distribution of various samples of close 
binaries with solar-type components confirmed the abovemen- 
tioned scenario, sugges ting as well a saturated braking law 
jRucinskill 19831 IMaceroni & Van 't Veerlll99ll IStepierilfl99l 
iMaceroni & Rucinskil 19991) . 

The expected physical properties of BW3 V38, as guessed 
in paper I solely on the basis of the photometry, put the system 
exactly at the proposed border of the efficient braking and make 
the system of paramount importance for the study of magnetic 
braking mechanisms and the effects of fast rotation in late-type 
dwarfs. Its location in the C-M diagram is also of great inter- 
est, as our system lies exactly midway between the two - and 
for a long time the only - other known eclipsing binaries with 



late dwarf components, YY Gem and CM Dra, fundamental 
lower-MS calibrators, and next to its recently discov ered, and 
presumably younger "twin", CU Cnc iDelfosse et all 1 999 i. 

Therefore, a spectroscopic study of this unique system was 
absolutely needed. In this paper we report the results of the ob- 
serving runs aimed to get mid-resolution spectra of the binary. 
Section 2 describes the details of data acquisition and reduc- 
tion, Section 3 presents the simultaneous solution of the light 
and radial velocity curves and a discussion of the system physi- 
cal properties, Section 4 gives an outline of the possible orbital 
evolution, followed by the conclusions. 



2. The spectroscopic observations 

BW3 V38 is a rather difficult target for radial velocity deter- 
mination. Its period requires short exposures to avoid Doppler 
smearing of spectral lines. Assuming, as a rule of thumb, that 
an acceptable exposure should last no more than 5% of the or- 
bital period (i.e. 0.05 in phase), the maximum duration shall 
not exceed 900 s . Therefore, because of its faintness, a medium- 
large class telescope is required. Furthermore, the field is rela- 
tively crowded and the binary has a nearby (at ~ 2" distance) 
and brighter (by ~ 2 mag ) optical companion. 

BW3 V38 was observed at ESO - La Silla during two half 
nights, in August 2001 (a previously scheduled run in July 2000 
was completely lost because of bad weather). A total of 17 long 
slit spectra were obtained, at various orbital phases, with the 
ESO 3.5m New Technology Telescope (NTT) and the ESO 
Multi-Mode Instrument (EMMI) in the red mid-resolution 
spectrographic mode. The slit width was fixed at 1". The spec- 
tra are 640 A wide and are centered around the Ha line (A = 
6563 A ). The selected grating (#6) has a nominal resolution of 
5500 (at 6000 A and with the l"-wide slit) and a dispersion of 
12.5 A mm" 1 (or 0.3 A pix" 1 ). The spectra of a few other K-M 
dwarfs, to be used as templat es, were acquired as well; the stars 
were c hose n from the lists oflKirkpa trick. H enry. & McCarthy! 
Jl99lh and lHenrv. Kirkpatrick. & Simons! Jl994l) . Finally, we 
selected the two templates that looked more similar to the spec- 
tra of the binary and of the optical companion, i.e. GJ 729 
(M3.5 V) and GJ 775 (K5 V), f or which accurate r adial ve- 
locities were recently obtained bv lNidever et alJ ll2002t) . 

Because of the above-mentioned constraints the exposures 
were limited to 900 s and the slit was aligned along the direction 
connecting our main target with the optical companion (and 
not - as usual - according to the parallactic angle). In this way 
the spectra of both objects could simultaneously be acquired 
and the contribution of the close-by star could be taken into 
account in data reduction, providing - as an extra bonus - a 
useful internal check. 

The data reduction was performed by means of standard 
techniques (the long slit reduction routines of the IRAF 1 pack- 
age). The steps were bias subtraction, flat fielding, (partial) 
cosmic ray removal, one-dimensional spectrum extraction and 



1 IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of the 
Universities for Research in Astronomy, Inc., under cooperative agree- 
ment with the National Science Foundation 
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Fig. 1. Three smoothed spectra of BW3 V38, from bottom to 
top: around primary eclipse, quadrature and secondary eclipse. 
The spectra are normalized to the highest continuum level. The 
splitting of Ha at quadrature indicates the presence of strong 
chromosphere in both components. 
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Fig. 2. Detail of the Ha region in the spectrum at phase (f> = 
0.233, note the partial reversal of the profiles. 

wavelength calibration. The radial velocities were obtained by 
standard cross correlation technique (as provided by the IDL 
5.5 and IRAF packages) after a few preliminary steps of: noise 
reduction (smoothing by convolution with a 3-pix FWHM 
gaussian curve, i.e. smaller than the slit projection), normaliza- 
tion to a pseudo-continuum (obtained by means of the IRAF 
continuum polynomial fitting and rejection algorithms) and 
transformation to the heliocentric system. The orbital phase de- 
termination s howed an excellent a greement with the ephemeris 
published by Udals ki et al.l dl995): in fact the same binary was 



briefly observed during the second phase of the OGLE exper- 
iment, a few years later, and no measurable period change or 
phase shift was found (Udalski, private communication). 

At this stage, an inspection of the radial velocities of the 
optical companion revealed the existence of a non-negligible 
nightly positive trend in the derived velocities (with a range 
of ^ 30 km s~'), corresponding to ~ 2 pixels. While its ori- 
gin remained unknown, but could presumably be ascribed to 
mechanical problems, the trend could be removed completely 
by a careful use of the night sky emission lines. The disper- 
sion solution was recomputed, using for wavelength reference 
the night-sky lines, which could be identifie d using the spec- 
tral at las of OH and O2 emission lines by lOsterbrock et alJ 
( 1996). Each spectrum was wavelength-calibrated with its own 
sky lines (with the exclusion of the bright template stars, whose 
exposure were too short to get a usable sky level; in that case 
the nearest exposure of the binary was used). In this way the 
trend completely disappeared and the standard deviation of the 
radial velocity of the optical companion decreased from 9.3 to 
2.2 km s _1 . 

Figure^shows three spectra of the variable close to the two 
minima and to quadrature. The typical S/N ratio of the binary 
spectra is ~ 10, i.e. still acceptable for radial velocity determi- 
nation. However, the task was not expected to be easy, because 
of the late spectral type and the related problems of continuum 
normalization, and because of the high rotational velocity of 
the components. 

On the basis of Figure ^ an d of the parameters we know 
from the light curve solution of paper I (in particular the high 
inclination that assures almost total eclipses), we can conclude 
that the binary is indeed formed by two active similar late-type 
dwarfs. While we do not have in our spectral range features that 
allow a quantitative spectral classificati on (as, for the red part of 
a late dwarf spectrum, the Ti05 index o flReid. Hawlev. & Gizisl 
Jl995l) 1. the general shape of the spectrum and the visibility of 
the TiO bands at 6569 and 6651 A suggest a M3 V classifica- 
tion. 

The second firm conclusion we can draw from the figure 
is that both components are dMe stars, as a splitting of Ha in 
emission is clearly visible at quadrature. For both components 
the Ha equivalent width around quadrature is 2.7 A, a 
value indicating the presence of strong chromospheres. Among 
dwarfs with Ha in emission, however, higher EW values are 
found, especially in component s of wider (and pr esumably 
younger) binary systems dGizis. Reid. & Hawlevll2002i) . typi- 
cal values span the range 0-10 A . 

The age of BW3 V38 is unknown, as well as its population 
characteristics; however its location in the galactic plane (pa- 
per I) and the low value of at least the barycentric velocity (see 
Table offer a weak support for a young/ intermediate disk 
population. Regardless, our system must be significantly older 
than the very similar but wider system CU Cnc (GJ 2069A). 
This too is formed by two M3 dwarfs, but has an orbital period 
of 2^77 and an estimated age, from the like ly associatio n with 
the Castor moving group, of only =; 300 My (Ribas 2003, here- 
after R03). According to lGizis. Reid. & HawlevN2002l) . the Ha 
EW of CU Cnc is = 4.6 A and its activity level, as estimated 
from the ratio L Ha /Lb i, is very high. R03 reports individual 
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Fig. 3. The radial velocity curve of BW3 V38. Dots indicate 
the primary and triangles the secondary velocity curve, the full 
line the result of the simultaneous solution, the dotted line the 
stellar barycentric radial velocities. Empty symbols correspond 
to points that were given half weight in the simultaneous solu- 
tion. The vertical bar is the mean uncertainty of the observed 
points. 
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Fig. 4. The new light curve solution BW3 V38 (lower panel) 
and the deviations. Paper I solution is virtually indistinguish- 
able from the new one and was not drawn, to avoid confusion. 





Fig. 5. The 3-D drawing at orbital phase 0.25 (upper panel) and 
the equatorial section (lower panel) of BW3 V38. The dashed 
line in the lower panel indicates the critical inner Roche lobe. 



were fitted with gaussian curves, leaving as free parameters the 
position and height. The amplitudes were estimated from the 
fit close to minimum phases (the light curve solution indicates 
almost total eclipses of two similar size components) and kept 
constant for the other orbital phases. As the least square fitting 
routine provided, as usual, unrealistically small errors, a value 
between 20 and 30 km s _1 (depending on the component and 
phas e) was estimated by means of the IRAF implementation of 
the lTonrv&Davisl ll979) algorithm, where the velocity errors 
are also computed, based on the fitted peak height and the anti- 
symmetric noise. This estimate tells us that the accuracy of the 
derived orbital elements will not be high enough to provide a - 
badly needed - improvement of the lower main sequence cal- 
ibration. However, a full exploitation of the existing data will 
improve our knowledge of this unique binary system. 



EW values for each component (of 3.85 and 4.05 A for the pri- 
mary and the secondary) and a value of LJLb n] = 10" 3 , that 
presumably marks the activity saturation limit ( Delfosseet al. 
1998). 

The shape of the Ha profile of BW3 V38, that around 
quadrature clearly shows a partial central reversal for both com- 
ponents (Figure indicates a lower activity level than typi- 
cally observed in most close binaries containing d Me dwarfs, 
which usually show a peaked profile ( Stauf fer & Hartm ann 
1986). 

Figure [3] shows the radial velocity curves, as derived by 
the cross correlation technique. The cross correlation functions 



3. Physical properties of BW3 V38 

A simultaneous solution of the radial velocity curves and the 
OGLE I-band light curve was computed by means of the 
Wilson and Devinney code ( I Wilson & Devinnevlll97lL in its 
2001 distribution 2 ). The input parameters for the light curve 
solution were taken from paper I and those for the radial ve- 
locity curves were obtained by a preliminary sine fit of the 
velocity curves. Looking for further information on the com- 
ponent color and temperature, we checked if the sys tem was 
detected in the 2MASS survey dSkrutskie et alJI 19951) . but the 

2 the code is kindly made available by the author at 
http://www.astro.ufl.edu/pub/wilson 
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Table 1. Simultaneous light and radial velocity curve solution 
ofBW3 V38 



Parameter 


this paper 


paper I 


i (degr) 


o c c 1 i r\ o a 

85.51 ± 0.84 


85.7 ± 1.2 


q 


0.950 ± 0.067 


0.77 ± 0.09 


£2, 


3.71 ±0.11 


3.48 ±0.13 


Q. 2 


4.02 ±0.19 


3.55 ± 0.25 


T 2 (T,=3500) 


3448 ± 11 


3459 ± 15 


L,/(Li+L 2 ) 


0.594 ± 0.031 


0.597 ± 0.036 


A<p 


0.0023 ± 0.0005 


0.0020 ± 0.0006 


tj (side) 


0.372 ± 0.008 


0.380 ±0.011 


r 2 (side) 


0.323 ± 0.028 


0.322 ± 0.050 


A(R ) 


1.355 ±0.066 




Ki (km s-») 


167 ± 8 




y (km s~') 


-0.6 ±5.9 





only source at the appropriate location corresponds to a K-type 
star, i.e. to the brighter optical companion (unfortunately its 
optical separation with the binary is comparable to the survey 
true PSF resolution of * 2.5" (FWHM)). We kept therefore the 
primary effective temperature derived in paper I. 

As in paper I, the photometric observations were given 
equal weights, while the radial velocities were unequally 
weighted (w, = 0.5 for phases close to the photometric minima, 
see Fig. [3}. The curve dependent weights (required to weigh 
the different sets of observations with respect to each other) 
were assigned, respectively, on the basis of the standard devi- 
ation of the preliminary sine fit for the radial velocity curves 
and of the local fit at maxima for the light curve. The corre- 
sponding values are cr„ = 25 km s _1 and cr\ c = 0.02 (in light 
units). The iterative solutio n used the multiple subset scheme 
JWilson & Biermannl 19761) . The adjusted parameters were the 
orbital inclination, i, the mass ratio 3 , q = m^lmx, the surface re- 
duced potentials, Q12, the secondary effective temperature^, 
the primary fractional luminosity, Li /(Lj + L2), the orbit semi- 
axis, A, the barycentric velocity y, and the phase shift A(p. The 
final solution is shown Figures 0} and the final parame- 
ter values are collected in Tabled The latter contains as well 
the component "side" fractional radii (i.e. the fractional radii 
in the orbital plane and at 90° from the axis connecting the 
star barycenters) and the primary radial velocity curve ampli- 
tude, Ki . The other non-adjusted parameters defining the so- 
lution maintained the same values of paper I, namely the lin- 
ear limb darkening coefficients x\£ = 0.56 at A = 7520 A 
(and X\ t 2 = 0.65 at Ha), the gravity darkening coefficients 
/?t,2 = 0.32, the albedos = 0.5. Moreover, the dark spot of 
paper-I was kept unchanged in the simultaneous solution and 
is defined by longitude and angular radius, respectively of 290 
and 1 1 degrees, and the temperature factor with respect to the 
unspotted surface, TJT* = 0.8. All the errors reported in the 
tables are standard errors. 

3 Star "1" is defined in the code as the component eclipsed at pri- 
mary minimum, in our case q < 1 and star "1", the primary, is as well 
the more massive component. 



Table 2. Derived physical quantities for BW3 V38 compo- 
nents. 



Parameter 


primary 


secondary 


(solar units) 






in 


0.44 ± 0.07 


0.41 ± 0.09 


R 


0.51 ±0.04 


0.44 ± 0.06 


M bol 


8.39 ± 0.08 


8.78 ±0.13 


logg 


4.66 ±0.13 


4.77 ± 0.21 



based on M bo] = 4.75 



The introduction of the radial velocity curves in the solution 
slightly changes the system configuration (Fig. [3 with respect 
to that obtained in paper I. The most affected parameter is - as 
one could expect - the mass ratio, that increases from q = 0.77 
to q — 0.95. The other adjusted parameters vary as predicted by 
the correlation plots of paper I, that were computed by means 
of bootstrap "resampling". The configuration is still a detached 
one, but the primary is slightly closer to the critical lobe (the 
ratio between the mean stellar radius and that of the critical 
surface increases from r/r« = 0.95 to r/r« = 0.98), while the 
secondary is somewhat more inside ( r/>« decreases from 0.90 
to 0.86). The increase in relative radius is balanced by the de- 
crease of secondary temperature, so that the fractional lumi- 
nosities are almost unchanged. 

It is well known that in detached and partially eclipsing sys- 
tems the light curve shape is only weakly dependent on the 
mass ratio and that there is a rather strong correlation between 
the geometric parameters. This explains why the maximum de- 
viation between the new and the paper I synthetic light curves, 
corresponding to the two solutions, is of only 0.005 mag (and 
the mean difference of just 0.001 mag). 

The derived absolute dimensions are presented in Table |2] 
The uncertainties on the physical parameters are unfortunately 
large, nevertheless the location of the components in the H -R 
and in the mass radius diagrams (Figure[6} are consistent with 
those of typical M-dwarfs. 

Figure [6] shows a comparison between the location of the 
known double-lined eclipsing binaries with M dwarf com- 
pon ents and three theoretical isochrones of different a uthors 
(lD , Antonalll2004 : lBaraffe etaD dl998t) : IVi et alJll200ll) ). 

The binary s ystems are: the already mentione d CU Cnc 
(R03), CM Dr a kacvlll 977t iMetcalfe et alJIl 996h YY Gem 
and V8 1 8 Tau JTorres & Ribasll2002h . whose secondary com- 
ponent is an M dwarf. It has to be stressed that CM Dra is a 
non-homogeneous system with respect to the rest of the sam- 
ple (and to the th eoretical iso chrones), as it seems to be a pop- 
ulation II object llLacvll 19771) . The physical parameters of the 
various objects are taken from the mentioned papers. 

The isochrones are for an age of 1 Gy and solar chemical 
composition; our system is probably older, but evolution affects 
only slightly the isochrone location for the mass range of inter- 
est (the isochrone age was actually chosen as a compromise 
between the age of the other younger systems and our object). 
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Fig. 6. The location of BW3 V38 components (points with error bars) in the mass radius and in the logT e - My diagrams. 
The locations of the other known eclipsing binaries with M-dwarf component(s) are also shown: the filled dots indicate CU Cnc 
components (R03), the starre d dot the secondary star of V818 Tau, the square YY Gem mean conponent (TR02), the triangles 
CM Dra components JLacvll977llMetcarfe et aljl 996). The theoretical isochrone s (1 Gy, solar c hemical composition) are from 
lD , Antonald2004 (continuous line). lBaraffe et alJ ( 1 1 9981) (dashed-dotted line) and lYi et alJd200ll) (dashed line). 



All the other systems but CM Dra are believed indeed to be 
younger: CU Cnc because of the association, on the basis of 
its space motion, with the Castor moving group, to which YY 
Gem presumably belongs, and V8 1 8 Tau as being a member of 
the Hyades cluster. 

Though the physical parameters of BW3 V38 are not ac- 
curate enough to pose constraints on the lower main sequence 
calibration, its location on the mass-radius diagram is - gener- 
ally - consistent with the conclusions drawn from all the other 
systems, i.e. that the theoretical models seem to underestimate 
the MS stellar radii and/or overestimate the effective temper- 
atures. A comparis on with a larger set of models, performed 
by Ribas (R03) and Torres & Ribas (2002) shows that this un- 
derestimate can be as high as 20 % for the stellar radius, and 
that most models overestimate temperatures (by 150 K). This 
is regardless of the recent improvements in the input physics 
for low mass stars. (The apparent agreement, in Fig [6] of CM 
Dra parameters is obviously due to the difference in chemical 
composition of the object with respect to that of the theoretical 
isochrones; lower metallicity isochrones would correspond to 
smaller radii). 

The radiative properties of BW3 V38 seem in good agree- 
ment with those of the other dwarfs in the same mass range, 
however, at variance with CU Cnc, the temperatures could not 
be independently derived. 

Both BW3 V38 components are very fast rotators: the radii 
of Table [2] and the hypothesis of spin-orbit synchronization 
yield equatorial velocities of, respectively, v\ = 131 + 10 and 
v 2 = 113 ±15 km s" 1 . 

iMullan & MacDonaldl feOOll) propose that the larger radii 
and lower temperatures of observed M dwarfs are due to 
the presence of strong magnetic fields. There is the indi- 
cation, indeed, that active dMe dwarfs (V-I < 2.7) are 
brighter on average by about 0.5 mag than dM stars of the 



same color faawlev. Gizis. & Reidl [l996), and that the ex- 
cess is even larger in clusters: 0.7 mag in the Hyades and 
1.2 mag in IC 2602 . The non-standard stellar models of 
Mul lan & MacDonaldl feOOll) . that include a simplified treat- 
ment of magnetic field effects, agree with these observations, 
suggesting larger radii and lower temperatures for stars host- 
ing strong fields (up to several tens of MGauss at the base of 
the convective zone). In addition in their models the limit of 
transition to a fully convective structure is shifted to masses 
lower than the standard model typical value 0.35/0.40 M Q . As 
a consequence the aQ. dynamo can still be at work for these 
mass values, and the AML through magnetic braking as well. 
The reason of the dearth of close systems with M type com- 
ponents shall therefore be found in the timescale of the AML 
process, rather than in its inhibition. 



4. An outline of BW3 V38 orbital evolution 

With the physical parameters of Table|2]the angular momentum 
content of BW3 V38 is H oth = 1.3 • 10 51 in cgs units (without 
the spin AM that amounts anyway to < 10 5l) cgs, i.e. less than 
10%). This is presumably only a fraction of the original con- 
tent. 

There are good reasons to believe that systems like 
BW3 V38 form with a much larger separation and evolve, later 
on, toward closer orbit s. The current knowledge of close bi- 
nary formation (Ghez 1996; B onnellll200ll) rules out the fis- 
sion process that does not seem able to produce stable binary 
configurations. The other "in situ" formation process, capture 
and hardening of the binary by ejection of a component from 
a triple system, is rare and inefficient in the field and works 
against the presence of small mass stars in the close pair, as it 
is generally the less massive star that is ejected. 



Table 3. Orbital evolution parameters 
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1.6 



Model 


a 


k 2 


i 


tm (Gy) 


A 


0.25 


0.2 


30° 


2.7 


B 


0.25 


0.2 


90 3 


4.5 


C 


0.25 


0.1 


90° 


8.9 


D 


0.00 


0.1 


90° 


9.8 



The hydrodynamical calculations of binary formation by 
fragmentation of Bate, Bonnell & Bromm (2002) indicate that 
proto-binaries form with large separations (> 10 AU) and go 
through a phase of accretion and orbital evolution towards 
tighter orbits ( ~1 AU). Subsequent less known phases of or- 
bital evolution, presumably due to interaction with the matter 
still present in the system, should finally produce really close 
binaries. 

On observational grounds, lMeloetalJd200ll) report in their 
collection of data for known PMS spectroscopic binaries that 
the shortest orbital period (relevant to a G5+K1 system) is P = 
1<?67. 

As a thorough description of the formation processes from 
the initial cloud to the Pre Main Sequence (PMS) is not yet 
available, a rough estimate of a minimum separation can be de- 
rived by the simple argument of the orbital semi-axis necessary 
to accommodate the pre-main-sequence (PMS) progeni tors of 
the binary. The PMS tracks of iD'Antona & Mazziteliil dl994 
for a 0.4 M star give a radius, during the D-burning phase, of 
~ 3 R Q , corresponding to a period of ~ 1^5 . 

Several authors have computed the orbital period evolution 
of a synchronized binary under the influence of AM L (Vilhu 
1982 [lRucinskll983t|Ma"ceroni & Van 't Veerl 199 it iMaceronil 
19921: IStepieiil 119951 llvlaceroni & Rucinskilll999l) developing 
an initial suggestion bv lVan't Veerl 111 97911 . Most works focused 
on deriving the slope of the braking law, without directly com- 
puting the amount of AML. Moreover, the greatest part of the 
abovementioned works concerns the better known case of sys- 
tems with solar-type components. We can try at least to es- 
timate the expected evolution of the orbital period using the 
available data on M dwarfs. 

The spin-down of a single late-type star can be described 

by: 



"Csd OJ 



(1) 



where c sd and a are constants and u> is the angular velocity. An 
elementary integration yields (a + 1): 



OJ = oj 1 



t-t 

(1 -a)T sd 
where the quantity 



l/Cl-a) 



Csd 



(2) 



(3) 



is the spin-down timescale (for t = f„). In a close binary thanks 
to synchronization (which can be considered as instantaneous, 




age (Gy) 

Fig. 7. The initial period, as function of system age, that is re- 
quired to reach by AML the present period of BW3 V38, ac- 
cording to the scheme proposed in the text. The curves corre- 
spond to different values of the parameters entering the compu- 
tation, as indicated in Table [3] 

taking place on a much shorter timesc ale than AML, see for in- 
stance M aceroni & Van 't Veerl (1 1991)) the loss of spin angular 
momentum from both components: 



k 2 {mir\ + rri2r\) to 



-k (m\r, + m 2 r 2 ) c sd of 



can be equated to the orbital angular momentum loss: 
1 



orb 



,2/3 



/i , \-2 5/3 -4/3 ■ 

q ( 1 + q) mJ a> ' oj 



(4) 



(5) 



with mi being the total mass. It is assumed that <w 01 -b = co sp = co, 
and that the masses, the radii and the gyration radius, k 2 , do not 
appreciably change with time. Besides, the small contribution 
of the spin to the total AM (the quantity to which Equation [5] 
should apply) is neglected. 

From Eqs d?]|5}, and introducing the orbital period, we get: 



P = - 
with 



c P 



,2/3-a 



3 (27r) ff+1/3 k 2 ( mi r\ + m 2 r 2 2 ) c sd (1 + q) 2 
G 2/3 ml! 3 q 



(6) 



(7) 



Equation l|6) can be integrated once suitable values for a, c sd 
an d k 2 are chosen 

iDelfosseet all i ll 998b derive that the spindown timescale of 
field M dwarfs at spectral type M4 is "a significant fraction of 
the age of the young disk (~3 Gy)", a reasonable choice for M3 
components is T sd 1 Gy. 

To derive c sd from Eq.|5]we also need a value for a> . The 
abovementioned paper gives the observed shortest rotation pe- 
riod {PI sin i) of a field M4 dwarf belonging to the young disk 
population (G 165-08), its value being about 7 h . This value is 
derived from vsirn measurements, but it is reasonable to as- 
sume that the fastest observed rotation corresponds as well to 
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high inclination. In the following we also computed, as a check, 
an extreme case with a very low inclination value. 

To choose a we recall that both the studies on single 
stars and on binaries with late type components indicate sat- 
uration of the braking law. The saturation boundary at v > 
10 km s _1 corresponds, with the radii of BW3 V38 compo- 
nents, to an orbital period ~ 2^6, meaning that most if not 
all of its peri od evolution must have take n place in a satu- 
rated regime. Ma ceroni & Rucinskil dl999h . from a study of 
the period distribution of late type binaries with equal compo- 
nents selected from the OGLE I databa se, find a "best value" at 
a = 0.25. Previous studies, (as Maceroni & Van 't VeeJ Jl99ll) 
and lMaceronil i 1 9921) '). suggest a 0. The models and the anal- 
ysis of M dwarf data suggest similar slopes, with, this in case, a 
shift of th e saturation boundary to long er periods with decreas- 
ing m ass iSills. Pinsonneault & Terndrupl200ollPizzolato et alJ 
l2003h . Finally, the value of k 2 does not vary much with mass 
along the main sequence, going from k 2 = 0.1 for a so- 
lar type star to k 2 = 0.2 f or a fully convective brown dwarf 
jArmitage & BonnelfeOO^ . 

Equation|5]was integrated with the different sets of parame- 
ters of Tableland the results are shown in Figure^ The curves 
were obtained integrating the equation backwards in time and 
taking as a boundary condition the present period of BW3 V38. 
The figure gives, therefore, as function of the system age, the 
initial period from which our system could have evolved to the 
present configuration. The extreme case with i = 30° corre- 
sponds of course to the fastest evolution. The value of r s d was 
kept fixed, but any larger value would make the evolution still 
slower. 

In all cases period decrease is rather slow. The evolution of 
a system similar to CU Cnc (i.e., P m = 2 ( ?77) would take, to 
reach the BW3 V38 configuration, a time larger than the age of 
the Universe (assumed to be 14 Gy). 

One may wonder if other AML mechanisms, such as AML 
by gravitational radiation, could have significantly contributed 
to the orbital evolution of the system. The time-scale of pe- 
riod change due to grav itational waves can be estimated from 
JPress & Thome Il972l) : 



2.8 x 10 



7 (mi + mj) 
m\ni2 



1/3 



3 8/3 



(8) 



with tq in years and P in hours. For an initial period of the 
order of a day (and using the masses from Table [2), we get 
tg - 7x10" y, while, for the current period value, tg - 9.4 Gy. 

For comparison purposes, a timescale of AML via magnetic 
braking, tm = \P/P\, can readily be derived from Eq.[6]and its 
value, for the current period, varies in the range 2.7 - 9.8 Gy 
(last column of Table[3j. The contribution of gravitational radi- 
ation could, therefore, be relevant in the future orbital evolution 
of the system, and only in the cases where the magnetic braking 
has the minimum efficiency. However - given the steep depen- 
dence of tq on P - its contribution was negligible during the 
past evolution of the binary. 

In conclusion, in the outlined scenario, if BW3 V38 has an 
age typical of an intermediate disk population (say, 5 Gy) its 
initial period could not have exceeded P m ~ l d and, even with 



all the uncertainties related to the AML description, it is not 
surprising that systems similar to it are very rare. 

5. Conclusions 

The analysis of the spectroscopic data of BW3 V38, in spite of 
their moderate quality, significantly increases our knowledge 
of this interesting system. 

The classification of the components as M3 dwarfs (which 
was previously based only on the V-I color from OGLE-I) is 
now confirmed. Both components are indeed dMe stars, with 
strong chromospheres, as indicated by the Ha emission pro- 
files. The Ha equivalent widths suggest an average activity 
level, when compared to similar binaries. This is an indica- 
tion that the system is older than similar systems. Thanks to 
BW3 V38 we have clear evidence that binaries with M dwarfs 
and very short periods are rare objects but do exist. If they form 
by AML from wider orbits, the AML timescale will then be < 
the age of the Galaxy. 

The simultaneous radial and light curve solution provided 
a consistent set of absolute elements. The components are very 
similar stars, with the slightly more massive one very close to 
the limiting Roche surface. The system can be considered a 
tighter twin of CU Cnc. 

The derived masses and radii agree - within the limits of the 
large uncertainty - with the relations obtained from the M-type 
components of the other known systems, which indicate a dis- 
crepancy, of the order of 10%, between the theoretical models 
and the values determined from eclipsing binary stars. It has 
been suggested that the effect could be related to the presence 
of strong magnetic fields, as dMe stars are indeed brighter than 
the inactive dM dwarfs. In addition the BW3 V38 components 
are very fast rotators (130 km s _1 ) and the larger radii could 
also be related to this property. 

A simple estimate of the orbital evolution of the system, 
taking into account the latest findings on the spin down of late 
type dwarfs, shows that the AML loss process is rather slow 
(the evolution from a configuration similar to that of CU Cnc 
would take, in the fastest case, ~ 14 Gy). In this frame the 
rareness of binaries similar to BW3 V38 can easily work. 

It is clear that the system deserves better quality observa- 
tions. It could be a fundamental system to better constrain the 
lower MS calibration, and to study the effects of magnetic fields 
and fast rotation, but higher S/N spectra are needed. These will 
require a telescope of the 8m class. We hope that the results 
obtained with this study will trigger the interest in this unique 
object, and qualify it as a priority target for the immediate fu- 
ture. 

Acknowledgements. The collection of data for this system was par- 
ticularly difficult, as weather and technical problems impaired several 
efforts in this direction. CM especially thanks Slavek Rucinski, for 
continuous encouragement and for his help during a first attempt at 
CFHT, that unfortunately did not produce any useful photons. The 
team of ESO-NTT is also acknowledged for precious advice during 
and after the two observing runs at ESO. The authors are also in- 
debted to Vincenzo Testa and Gian Luca Israel for discussions on 
data reduction techniques, to Francesca D'Antona for providing new 
isochrones and to the anonymous referee for useful comments. This 



C. Maceroni and J. Montalban: Absolute elements of BW3 V38 



9 



research project was funded by MIUR/Cofin and F-INAF programs 
and by ESA-Prodex under contract 1 5448/0 l/NL/SFe(IC)-C90 135. 



References 

Alvarez, R. & Plez, B. 1998, A&A, 330,1109 
Armitage, P. J. & Bonnell, I. A. 2002, MNRAS, 330, LI 1 
Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, 
A&A, 337, 403 

Bate, M. R., Bonnell, I. A., & Bromm, V. 2002, MNRAS, 336, 
705 

Bonnell, I. A. 2001, in The Formation of Binary Stars, IAU 
Symp. 200, ed. H. Zinnecker & R. D. Mathieu (Kluwer), 23. 

D'Antona, F. & Mazzitelli, I. 1994, ApJS, 90, 467 

D' Antona F. 2004, private communication. 

Delfosse, X., Forveille, T., Perrier, C, & Mayor, M. 1998, 
A&A, 331,581 

Delfosse, X., Forveille, T., Mayor, M., Burnet, M., & Perrier, 

C. 1999, A&A, 341, L63 
Durney, B. R., De Young, D. S., & Roxburgh, I. W. 1993, 

Sol. Phys., 145, 207 
Ghez, A. M. 1996, in NATO ASIC Proc. 477, Evolutionary 

Processes in Binary Stars, 1 
Gizis, J. E., Reid, I. N., & Hawley, S. L. 2002, AJ, 123, 3356 
Hawley, S. L., Gizis, J. E., & Reid, I. N. 1996, AJ, 1 12, 2799 
Henry, T. J., Kirkpatrick, J. D., & Simons, D. A. 1994, AJ, 108, 

1437 

Kim, Y. & Demarque, P. 1996, ApJ, 457, 340 
Kirkpatrick, J. D., Henry, T. J., & McCarthy, D. W. 1991, ApJS, 
77,417 

Lacy, C. H. 1977, ApJ, 218, 444 

Maceroni C. 1992 in ASP Conf. 40, Inside the star, eds W.W. 

Weiss and A. Baglin, 374 
Maceroni, C.,& Rucinski, S. M. 1997, PASP, 109, 782 
Maceroni, C. & Rucinski, S. M. 1999, AJ, 118, 1819 
Maceroni, C, Van 't Veer, F. 1991, A&A, 246, 91 
Maceroni, C, Van 't Veer, F. 1993, A&A, 277, 515 
Melo, C. H. F, Covino, E., Alcala, J. M., & Torres, G. 2001, 

A&A, 378, 898 

Metcalfe, T. S., Mathieu, R. D., Latham, D. W., & Torres, G. 

1996, ApJ, 456, 356 
Mochnacki, S. W. 1981, ApJ, 245, 650 
Mullan, D. J. & MacDonald, J. 2001, ApJ, 559, 353 
Nidever, D. L., Marcy, G. W., Butler, R. P., Fischer, D. A., & 

Vogt, S. S. 2002, ApJS, 141, 503 
Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., 

& Vaughan, A. H. 1984, ApJ, 279, 763 
Press, W. H. & Thorne, K. S. 1972, ARA&A, 10, 335 
Osterbrock, D. E., Fulbright, J. P., Martel, A. R., Keane, M. J., 

Trager, S. C, & Basri, G. 1996, PASP, 108, 277 
Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., & Ventura, 

P. 2003, A&A, 397, 147 
Reid, I. N., Hawley, S. L., & Gizis, J. E. 1995, AJ, 110, 1838 
Ribas, I. 2003, A&A, 398, 239 
Rucinski, S. M. 1983, The Observatory, 103, 280 
Rucinski, S. M. 1992, AJ, 103, 960 
Rucinski, S. M. 1997, AJ, 113, 407 



Sills, A., Pinsonneault, M. H, & Terndrup, D. M. 2000, ApJ, 
534, 335 

Skrutskie, M. F. et al. 1995, Bulletin of the American 

Astronomical Society, 27, 1392 
Stauffer, J. R. & Hartmann, L. W. 1986, ApJS, 61, 531 
Stepieh, K. 1995, MNRAS, 274, 1019 
Torres, G. & Ribas, I. 2002, ApJ, 567, 1 140 
Terndrup, D. M., Stauffer, J. R., Pinsonneault, M. H., Sills, 

A., Yuan, Y, Jones, B. F, Fischer, D., & Krishnamurthi, A. 

2000, AJ, 119, 1303 
Tonry, J. & Davis, M. 1979, AJ, 84, 151 1 
Udalski, A., Szymahski, M., Kaluzny, J., Kubiak, M., Mateo, 

M. & Krzeminski, W. 1995, AcA, 45, 1 
Van't Veer, F. 1979, A&A, 80, 287 
Vilhu, O. 1982, A&A 109, 17 
Wilson, R. E., Devinney, E. J. 1971, ApJ, 166, 615 
Wilson, R.E., Biermann, P. 1976, A&A, 48, 349 
Yi, S., Demarque, P., Kim, Y, Lee, Y, Ree, C. H., Lejeune, T, 

& Barnes, S. 2001, ApJS, 136, 417 



